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Abstract 

We consider Golonib rulers and their construction. Common rulers feature marks 
at every unit measure, distances can often be measured with numerous pairs 
of marks. On Golomb rulers, for every distance there are at most two marks 
measuring it. The construction of optimal — with respect to shortest length 
for given number of marks or maximum number of marks for given length — is 
nontrivial, various problems regarding this are NP-complete. We give a simplified 
hardness proof for one of them. We use a hypergraph characterization of rulers 
and Golomb rulers to illuminate structural properties. This gives rise to a 
problem kernel in a fixed-parameter approach to a construction problem. We 
also take a short look at the practical implications of these considerations. 
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Chapter 1 

Introduction 



A Golomb ruler is a specific type of ruler: Whereas common rulers have marks 
at every unit measure, a Golomb ruler only has marks at a subset of them. 
Precisely, the distance measured by any two marks on a Golomb ruler is unique 



on it. An example can be seen in Figure 1.1 Golomb rulers are named after 



Professor Solomon Golomb. According to various sources [8l[T2], he was one of 
the first to study their construction. 

Golomb rulers have various applications ranging from radio astronomy to 
cryptography. This explains the interest in computing Golomb rulers that are 
particularly short for a given number of marks or have many marks when given 
a maximum length. Unfortunately, from a computational complexity point 
of view, some of the corresponding decision problems have been proven to be 
NP-complete, while little is known about other very natural problems. 

Despite this, much effort has been made to compute short or dense Golomb 
rulers and to prove them optimal. Various implementations of exhaustive searches 
have been given and discussed as well as heuristic and evolutionary approaches. A 
sophisticated project searches for Golomb rulers through a distributed computer 
network, enabling users to donate idle computing time. 

In this work, we give a short insight into the work that has been done in the 
field and we briefly consider two natural problems of unsettled computational 
complexity. We give a natural hypergraph characterization for rulers such that 
only Golomb rulers correspond to a specific subset of the graphs. We then 
consider a construction problem that has been proven to be NP-complete. We 
give a simplified proof for this and then look at two natural parameterizations. 
For one of the parameterizations, we provide a fixed-parameter algorithm, and 
some heuristic improvements along with a cubic-size problem kernel that mainly 
follows from some structure that we observe in characteristic hypergraphs. 
Finally, we implemented an algorithm that uses the fixed-parameter approach 
and comment on our experimental results. 



1.1 Origins and Applications 

According to Colannino 8J and Dimitromanolakis [12 , W. C. Babcock first 
discovered Golomb rulers while analyzing positioning of radio channels in the 
frequency spectrum. He investigated inter-modulation distortion appearing in 
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Figure 1.1: A common ruler (left) and a Golomb ruler (right) is shown. For 
every distance in Golomb rulers, there is at most one pair of marks that measure 
this distance. For example, the distance one is only measured by the marks 
and 1 on the Golomb ruler whereas this distance is measured by six pairs of 
marks on the common ruler. Both rulers measure every integer distance up to 
their length. Such rulers are called perfect. This, however, is a rare trait among 
Golomb rulers, as one can proof |12) . 



consecutive radio bands [3] and observed that when positioning each pair of 
channels at a distinct distance, then third order distortion was eliminated and 
fifth order distortion was lessened. 

Rankin [29] lists other interesting applications, two of which we touch shortly 
here. In radio astronomy, arrays of radio telescopes are used to gather information 
about celestial bodies via interferometry. The telescopes are arranged in a single 
line, and information is extracted from difference measurements between two 
telescopes [S] . By placing them at the marks of a Golomb ruler, the number 
of these measurements and thus information gathered is maximized. This is a 
special case of a linear array. Linear arrays are also used in other related fields 
such as antennae construction. 

In computer communication networks Golomb rulers can be used to simplify 
the message passing process. When allocating the node names corresponding 
to marks on a Golomb ruler, messages do not need to specify both origin and 
destination addresses. Since the differences between marks in a Golomb ruler are 
unique, the difference and the direction of arrival suffice to identify the origin 
and destination node [5]. 



1.2 Preliminaries and Problem Definitions 

We now gather a common knowledge base for our considerations. At first we will 
introduce rulers and Golomb rulers, then go on to hypergraphs, some basic fixed 
parameter techniques and finally define some notation. We assume the reader 
to be familiar with basic mathematics and classic complexity theory. There are 
many recommendable books on complexity theory, see for example [21 [57] . 



1.2.1 Rulers 

Definition 1.2.1 (Ruler). A ruler is a set R := {nii : 1 < i < n} C Z 
with rrii < rrii^i. The mark rrii is called the i'th mark on R. The ruler is said 
to have n marks and |?7i„ — mi\ is called the length of the ruler. We call a 
set R' (- R a subruler of R. 

Definition 1.2.2 (Golomb ruler). A ruler R = {mj : 1 < i < n} is called 
Golomb ruler if for every (i G N \ {0} there is at most one solution to the 



equation d = rui — mj , mi 



eR. 



It is easy to see that if a ruler {rrii : 1 < z < n} is Golomb, so are the 
rulers {rriiC : 1 < « < n} and {rrii + c : 1 < i < n} for a constant c G Z. 

So for every Golomb ruler R = {rrii : 1 < * < n}, there is a Golomb ruler R' 
with only positive marks, starting with the mark and having the same set of 
differences {rrii — rrij : mi 7^ rrij A mi,mj e R}. A ruler R' is also called the 
canonical form of R |12j . 

Definition 1.2.3 (Perfect ruler). A ruler R = {mi : 1 < i < n} is called perfect 
if for every integer 1 < d < {rUn — mi) there is at least one solution to the 
equation d — mi — mj, rrii, ruj G R. 

It is not hard to see that there are no perfect Golomb rulers with more than 
four marks [29]. In fact, every Golomb ruler with n > A marks has length greater 
than n{n — l)/2. 

This insufficiency leads to following problem: 

Optimal Golomb Ruler 

Instance: n e N. 

Task: Find a Golomb ruler with n marks and minimum length. 

A natural decision problem (as posed by Meyer and Papakonstantinou [5^) is: 

Golomb Ruler Decision 

Input: n, Z? e N. 

Question: Is there a Golomb ruler with at least n marks and length at 

most Dl 



Definition 1.2.4 (Optimal Golomb ruler and the function G{n)). A Golomb 
ruler with n marks is called optimal if it is of shortest possible length. For every 
n £ N the function G{n) is defined as the length of an optimal Golomb ruler 
with n marks. 

No closed form expression is known for G{n) [12], and even using massive 
computational power, G{n) to date is only known for n < 26 [14] . 

1.2.2 Graphs and Hypergraphs 

A hypergraph iJ is a tuple (V, E) , where V^ is a finite set and E \s a, family of 
sets e such that e C 2^ \ {0}, where 2^ denotes the powerset of V . The elements 
of V are called vertices and the elements of E edges. Two vertices are called 
adjacent^ if there is an edge that contains both of them. A vertex v £V and an 
edge e € E are called incident, if w G e. 

Two hypergraphs H = {V,E) and H' — {V',E') are called isomorphic and 
we write H = H\ if there is a bijection (f) : V —> V such that the following holds: 

{wi,...,wje E^{<l){vi),...,(j}{vi)} eE' 

The hypergraph H' = {V' , E') is called subhypergraph or short subgraph of a 
hypergraph H — {V,E), ii V (^ V and E' <Z E. The hypergraph H is then also 
called a supergraph of H' . The hypergraph -ff [1^'] := {V, E') is called (vertex) 
induced subgraph of a hypergraph H = {V, E), ii V' ^ V and E' contains every 



set e, such that e € E and e C V' . A hypergraph Mh is called a minor of a 
hypergraph 7?, if Mh can be obtained from H by removing vertices, removing 
edges and contracting edges. Contracting an edge e means to delete e and every 
vertex contained in e and to introduce a new vertex v that is added to every 
edge that was incident to a vertex in e. 

The hypergraph H is called simple, if e 7^ e for every e, e G ii^. The 
hypergraph iJ is called d-uniform, for an integer d, if |e| = d for every e £ E. 

We draw a hypergraph by drawing points for every vertex, drawing curves 
for every edge e, |e| > 2, encircling all vertices in that edge, drawing straight 
lines between the contained vertices for edges e, |e| —2 and drawing loops at the 
contained vertex for edges e, |e| — 1. 

A graph is a simple 2-uniform hypergraph. A graph is called complete or a 
Clique, if all possible edges are present in the graph. A graph G = (Vi U V2, -E) is 
called bipartite, if Vi n V2 = and there is no e G i<^ such that e C Vi or e C V2 . 
A bipartite graph G — (Vi UV2,E) is called complete or a .ft'iViUVal i^ i^ contains 
all possible edges and still maintains its bipartite property. A graph is called 
planar, if it can be embedded in the plane, i.e. it can be drawn in a plane such 
that the drawings of the edges intersect only in their endpoints. A well known 
theorem by Kuratowski states that a graph is planar, if and only if it does not 
contain a clique with five vertices or a K33 as minor [34] . 

The incidence graph of a hypergraph H — {V, E) is defined as the bipartite 
graph I ^ {V\JE, E') with E' = {{v, e} : v e V Ae e E Av e e}. A hypergraph 
is called planar, if its incidence graph is planar. 

1.2.3 Basic Fixed-Parameter Techniques and Complexity 
Theory 

Many natural problems are NP-hard and thus are believed not to be solvable 
within running time bounded by a polynomial function. However, in practise the 
phenomenon can be observed that some instances of NP-hard problems in fact 
can be solved within reasonable time. This is because in classic computational 
complexity mostly worst case running times depending only on the input length 
are contemplated. The exponential worst case running times notwithstanding it 
is possible to identify structures that can be exploited by algorithms in some 
problems. Fixed-parameter algorithmics can be seen as the approach to find 
efficient algorithms not only with respect to input length but also to such 
structures — called parameters. We only recapitulate some very basic definitions 
and techniques here, for more on the topic see e.g. [211 UHl HHj- 

A parameterized problem is a language L C E* x S*, where S is a finite 
alphabet. The second component is called the parameter of the problem. In our 
problems the parameter will always be a nonnegative integer and therefore we 
restrict this definition to languages L C E* x N in this paper. A parameterized 
problem L is fixed-parameter tractable with respect to the parameter k if there 
exists an algorithm that decides L in f{k)p{n) time, where / : N — >■ N is a 
computable function only depending on k and p is a polynomial. 

Bounded search trees are a standard way to classify a problem as fixed- 
parameter tractable. Search trees are a way of systematic exhaustive search. 
The strategy is to find a small part of the input in polynomial time, such that 
at least one element of this part has to be in the solution. Then we branch 
into all cases of choosing one element of this part and recurse until a solution 



is found. The graph with the (recursive) caUs of the algorithm as nodes and 
an edge between two nodes, if one was called by the other is called the search 
tree. If in every node time is spent that is bounded by a polynomial in the 
input length and if we can bound the number of succeeding recursive calls at 
one node and the height of the tree by a function that depends only on the 
parameter, we obtain a fixed-parameter algorithm. If the algorithm has an input 
of size s and branches into recursively solving instances of sizes s — di, ..., s — di, 
then (di, ...,di) is called the branching vector of this recursion. Since search tree 
algorithms often terminate early and are easily parallelized, this technique has 
applications in practise. 

Let L be a parameterized problem. A reduction rule is a mapping of in- 
stances (/, k) to instances (/', k') such that the following conditions hold: First, 
(/, k) ^ L 4^ (/', fc') e L, which is also called correctness of the rule. Second, 
length(/') < length(/) and k' < k. An instance is called reduced with respect to 
a reduction rule, if the rule cannot be applied to the instance anymore. That is, 
the image of the instance under the reduction function is the same as the instance 
itself. A reduction to a problem kernel is a reduction rule that can be computed 
in 0((length(/))'^) time for some constant c such that length(/') < g{k), where 
(7 : N — > N is computable and depends only on k. The function g is also called 
the size of the problem kernel. 

Let L, L' be two parameterized problems. A parameterized reduction from L 
to L' is a function r : S* x N — )■ E* x N, (/, fc) i-)- (/', k') such that r is computable 
in /(fc)p(length(J, k)) time, (/, k) & L -i^ (/', k') G L' and k' depends only on k. 
Here, / is a computable function depending only on k and p is a polynomial. 

A parameterized problem L belongs to the class W[t] if there is a parame- 
terized reduction from L to a weighted satisfiability problem for the family of 
circuits of weft — the maximum number of gates with unbounded fan-in on an 
input output path — at most t and depth at most some function of the param- 
eter k. For an introduction to the W-hierarchy see [18]. We only use the fact 
that W[l]-hard problems are assumed not to be fixed-parameter tractable. A 
parameterized problem can be shown to be W[l]-hard by giving a parameterized 
reduction from a W[l]-complete problem. For example, the following problem is 
W[l]-complete with respect to parameter k. 

Independent Set 

Input: A graph G = {V, E) and an integer k. 

Question: Is there a vertex subset S QV such that fc < jS*] and G[S] 

contains no edges? 



1.2.4 Conventions 

We denote the number of vertices of a (hyper-)graph by n and the number of 
edges by rra, where it is not ambiguous. Also, we denote an edge of size d in 
hypergraphs by d-edge. Hypergraphs whose edges are of size exactly three or four 
play a major role in our considerations. We call these graphs 3,4-hypergraphs. 

In proofs, we use the lightning-symbol i to indicate an exposed contradiction 
and thus a finished case of a reductio ad absurdum. 



1.3 Previous Work 

In this section, we give a quick overview of the research that has been done in 
the Golomb ruler field in the past. 

1.3.1 Golomb Ruler Construction and Bounds on G{n) 

Dimitromanolakis ^12, analyzed the Golomb ruler problem with respect to the 
older so called Sidon set problem. Informally, a Sidon set is a finite set 5 C N 
such that the sum of any two elements of S is distinct. It can easily be shown 
that this definition is equivalent to the definition of the Golomb rulers. Similarly 
to Golomb rulers, there is an optimization problem for Sidon sets that asks the 
following: Given n, what is the maximum cardinality of a Sidon set S C {1, ..., n}? 

This problem has been studied extensively. Lower bounds have been given 
by Lindstrom [19] (as claimed in |12]) and Erdos and Turan [17] via construc- 
tion strategies. These bounds then have been applied to Golomb rulers by 
Dimitromanolakis [T^], yielding the lower bound G{n) > v? — In^n + -y/n — 2. 

Other construction strategies for Sidon sets and thus Golomb rulers have 
been given indirectly. Singer [3T] discovered a method for generating a set of 
(7+ 1 residues modulo g^ + (7+ 1 that form a Golomb ruler, q being a prime power. 
Dimitromanolakis [12] and Ruzsa [30] claim that Bose [7] found a strategy to 
do the same for q residues modulo 5^ — 1 and Ruzsa [SOJ found one for p — 1 
integers such that their pairwise sums are all different modulo p{p — 1), p being 
a prime number. The former implies that G{q) < q^ — 1 for prime powers q. 

Still, for all integers the best known upper bound is G{n) < Ir? + n. This is 
due to a relatively simple construction described in [12] . 

Erdos conjectured an upper bound to be G{n) < n^ + c, with c G M being a 
constant. Dimitromanolakis [H] has computed relatively short Golomb rulers 
and thus showed with computer aid that G(n) < ri^ for n < 65,000. 

1.3.2 Computational Approaches 

Numerous efforts have been made to compute Golomb rulers and prove their 
optimality with computer aid. We can only cover some selected cases here. 

Rankin [53] has developed some exhaustive search algorithms, reviewed their 
empirical running time, and provided a parallel implementation. He has been 
able to compute and prove optimal Golomb rulers with 17 and 18 marks in 
1993. Using a cluster of one Sun SparcServer 1000 and twelve Sun SparcClassic 
workstations, it took about 840 hours to complete the run for the 18-mark ruler. 

Distributed.net [TJ] — founded in 1997 — attacks the Golomb ruler problem 
through a globally distributed computer system. A client is provided, which 
enables the user to donate idle computing time to the project. The participants 
were able to show that previously shortest known Golomb rulers of length 24 
through 26 were indeed optimal. With 124,387 participants during the run for 
the 25-mark Golomb ruler, it took 3,006 days to find the result. The run for the 
26-mark Golomb ruler allegedly took only 24 days. Unfortunately, the source 
code of the projects key elements are not publicly available because of security 
reasons. 

Golomb rulers have also raised some interest in the field of evolutionary 
algorithms and other heuristic techniques. For example, see articles by Tavares 



et al. [53], Cotta et al. [TT], Pereira et al. [55] and Cotta and Fernandez [lU] . 

1.3.3 Complexity Theory 

Surprisingly, given the number of implementations, seemingly little is known 
about the computational complexity of GoLOMB Ruler Decision and Optimal 
GOLOMB Ruler. See recent publications by Meyer and Papakonstantinou [22] 
and Ma and Yao [20] • 

Meyer and Papakonstantinou 22j have focussed on the construction of 
Golomb rulers and proved the following problems to be NP-complete: 

GOLOMB SUBRULER 

Input: A finite set S* C N and n eN. 

Question: Is there a Golomb ruler S' C S with at least n marks? 

Golomb Ruler Sum 

Input: A finite set T C N and D, n e N. 

Question: Are there elements ti, ..., U G T such that n — I < i and the 

ruler {J2i=i tj : I < k < i}U {0} is a Golomb ruler of length equal to D? 

Golomb Ruler Subset Distances 

Input: A finite set of interval lengths T C N and n E N. 

Question: Is there a Golomb ruler R — {mi : 1 < i < n] such that 

{\mj — mi\ : \ < i < j < n} CI Tl 

Ma and Yao ^ have reduced Optimal Golomb Ruler to a problem called 
Seed Optimization. Unfortunately, they also note that the complexity of this 
problem is unknown. 

However, there has also been some research on intuitively related problems, 
i.e.. Difference Cover and Turnpike: A set A c N is called a difference 
cover for a set F C N if for each y &Y there exist at least two elements a, 6 G A, 
such that y = a — b. Mereghetti and Palano [H] have proven that a polynomial 
time algorithm for the following problem would imply P = NP. 

Minimum Difference Cover 

Instance: A set y C N. 

Task: Find the minimum cardinality A C N, such that A is a difference 

cover for Y. 

The Turnpike problem is defined as follows: 

Turnpike 

Input: A multiset of n{n — l)/2 integer distances. 

Question: Is there a set of n points in N with the given distances? 

Meyer and Papakonstantinou [22] note that the complexity of Turnpike is 
unknown. Variants and special cases of it have been studied, some of which have 
been proven to be in P, some to be pseudo-polynomial time solvable and others 
to be NP-complete. See references in [22] . 



Chapter 2 

Algorithms and Complexity 



Meyer and Papakonstantinou [55] provided insight into the computational com- 
plexity of constructing Golomb rulers. However, there are still white spots in 
the map of complexity of problems related to Golomb rulers. In [Section 2.1| we 
consider one central white spot, and try to explain why it still has not been 
settled. 



In Section 2.2 we introduce the notion of characteristic hypergraphs for rulers. 
This new technique helps illustrating problems related to Golomb rulers and we 
also gather some structural insights based on these graphs. The characterization 
serves as base for our considerations in the succeeding sections: 

[Section 2.3| contains an alternative and simplified proof for the NP-hardness 
of Golomb Subruler. (The original proof has been given by Meyer and 
Papakonstantinou [22].) We also consider the fixed-parameter tractability of 
Golomb Subruler with two natural parameters in [Section 2.4| We give a 
positive result for one of them and prove a cubic-size problem kernel. 

2.1 Notes on the Complexity of Optimal Golomb 
Ruler 

Recall the definitions of Optimal Golomb Ruler and Golomb Ruler 
Decision: 

Optimal Golomb Ruler 

Instance: n e N. 

Task: Find a Golomb ruler with n marks and minimum length. 

Golomb Ruler Decision 

Input: n,D €N. 

Question: Is there a Golomb ruler with at least n marks and length at 

most D? 



Although some authors [121 ES] beheve that Optimal Golomb Ruler is 
computationally hard, to date there is no proven evidence. Even Meyer and 
Papakonstantinou [55^ , who focussed on the complexity of problems related to 
Golomb rulers, do not state a conjecture on whether Golomb Ruler Decision 



10 



is NP-hard or whether it is in NP. We now discuss why these questions seem to 
be difficult to answer. 

At first, observe that the encoding of the maximum length D and the 
minimum number of marks n has a heavy impact on the complexity of Golomb 
Ruler Decision. This is due to the following. 

Theorem 2.1.1 (Berman [4j). If there is an NP-complete language L and a 
polynomial function / : N — > N such that the following statement holds, then 
P^ NP. 

yien:\{xeL: length(a;) < l}\ < f{l) 

In Golomb Ruler Decision, let D and n be encoded with binary alphabet 
and let Unary Golomb Ruler Decision denote the same problem, but with 
the input encoded with unary alphabet. In every unary language, there are at 
most I + 1 words of length at most I. Therefore, the number of positive instances 
of length at most I in Unary Golomb Ruler Decision is clearly bounded by 
a polynomial function in I. Thus, by |Theorem 2.1?T| and under the assumption 
that P 7^ NP, Unary Golomb Ruler Decision cannot be NP-complete. 
The problem clearly lies in NP, because the trivial certificate — a Golomb ruler 
satisfying the conditions — is of length 0(rilogZ3), which is polynomial in n 
and -D, so it cannot be NP-hard. 

For Golomb Ruler Decision, that is, using binary encoding of the input, 
it is not even clear whether it is in NP. This is due to the fact that, if it is in NP, 
then there must be a certificate of polynomial size for every instance. As noted 
above, the trivial certificate is of length O(nlog-D), which now is exponential 
in the input length (logn -t- \ogD). Obviously, this does not imply that such 
a certificate cannot exist, but it seems intuitively plausible that it comprises a 
Golomb ruler or some notion of the positions of its marks. To encode this into a 
word of size polylogarithmic in D and n seems to be a difficult task. 

Because the number of marks n is the culprit to the exponentiality of the 
trivial certificate here, it might seem obvious to try and encode only n in unary 
and give D in binary. However, this makes [Theorem 2 . 1 . 1 1 applicable again: 

Observation 2.1.1. Let L he a language and let f he a polynomial-time com- 
putable function such that either f{w) —1- or f{w) ~ M , where M is the coding 
of a Turing machine that decides in polynomial time whether w Cz L or not. Fur- 
thermore, let Lnt — {w G L : /(w) =-L}. Then L and Lnt are polynomial-time 
equivalent. 

Proof. It is clear that the language Lnt can be decided in polynomial time with L, 
because L„t C L. 

To reduce L to L„(, one simply computes f{w) for w Cz L. If f{w) =_L then 
w G Lnt a-nd we can output w. Otherwise, simulate /(w) = M with input w and 
output Wno or Wyes if M does not or does accept, respectively, where Wno 4- ^ 
and Wyes G L (we can code a constant number of such words in the coding of 
the Turing machine computing this reduction). D 

The interpretation of [Observation 2.1.1| is that any problem L is polynomial- 
time inter-reducible with a problem Lnt which contains only "non-trivial in- 
stances" , that is, all instances of L except those that are known to be decidable in 
polynomial time and that can be classified to be so in polynomial time. Therefore, 
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Figure 2.1; Two rulers with the marks a,b,c and rf, respectively. To the left, 
we see that the marks a and b measure the same distance as c and d. We 
consider this to be a conflict with respect to Golomb rulers and model it as 
an edge {a, &, c, d} in the corresponding hypergraph. To the right we see a 
degenerated form of a conflict which leads to an edge with only three vertices. 



for the sake of complexity classification, we can exclude trivial instances in every 
problem. 

We know from Subsection 1.3.1 ('[Golomb Ruler Construction and Bounds 



on G{n) ' ) that the length of an optimal Golomb ruler with n marks is upper- 
bounded by i e O(n^). Thus, for a given number n of marks, there are only 
polynomially in n many lengths D that form instances which can not trivially 
be checked in polynomial time. Out of the words of length at most I in the 
language induced by Golomb Ruler Decision with n in unary, D in binary 
and without trivial instances, there can be at most I + 1 words that represent 
distinct values of n and thus, there can be at most t{l + 1) words of length at 
most I. This means that if n is given in unary in Golomb Ruler Decision 
and D is encoded binarily, a language defined by this problem that excludes 
trivial instances again satisfies the condition of [Theorem 2.1.1] 

Also, since many NP-complete problems have a notion of efficient self- 
reduction, an approach to get a hint on the complexity of Golomb Ruler 
Decision would be to search for such a procedure. Self-reduction is a procedure 
to compute an optimal solution for an optimization problem using an oracle for 
the corresponding decision problem. Unfortunately, Meyer and Papakonstanti- 
nou [22] note that it seems difficult to efficiently compute a valid Golomb ruler, 
given an oracle for GOLOMB Ruler Decision. 



2.2 Characterizing Golomb Rulers through Hy- 
pergraphs 

In this section we provide a simple hypergraph characterization of Golomb 
rulers and consider structural properties of the implied hypergraphs. The 
characterization serves as base for considerations in the succeeding sections. 



2.2.1 Hypergraph Construction 

We start with a ruler i? C N and construct a hypergraph with edges consisting of 
either three or four vertices. Every mark corresponds to one vertex in the graph 
and every edge to one conflict, i.e., to a distance that is measured by two pairs 



of marks (see Figure 2.1). Constructing such a graph can be done by simply 



iterating over every three and four-tuple of marks and checking whether the marks 
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in the tuple have equal distances. See also |Algorithni HypcrgraphConstruction} 



Algorithm HypergraphConstruction: Constructing a characteristic 
hypergraph for a given ruler 

Input: A ruler i? C N. 

Output: A hypergraph Hj^ = {R,E). 

Start with an empty hypergraph H; 

Let R be the set of vertices in H; 

for (a, b, c, d) £ R^, \{a, b, c, d}\ = 4 do 
|_ if |a — &| ~ \c — d\ then add the edge {a, b, c, d} to H 

for {a, b, c) £ R^,\{a,b,c}\ = 3 do 
|_ if |a — 6| = |6 — c\ then add the edge {a, b, c} to H 

return H: 



Definition 2.2.1. In the following sections, we denote the hypergraph con- 
structed from the ruler R by Hj^ — (i?, E) and call it characteristic hypergraph 
ofR. 

The construction of Hji is clearly computable in 0(|i?|^). 

Lemma 2.2.1. Let R be a ruler and Hn = [R,E) the corresponding character- 
istic hypergraph constructed from R using \Algorithm Hyper graphConstruction\ 
Then R is a Golomb ruler if and only if E — 0. 

Proof. We show that both directions of the following equivalent statement hold: 
The ruler R is not a Golomb ruler if and only ii E ^ 0. 

First assume that E is not empty. Thus, there is either an edge {a, b, c} or 
{a, 6, c, d} in i? and the corresponding equations |a — &| = |c— d| or |a — 6| = |&— c|, 
respectively, hold. In other words, R contains a number of marks that have 
pairwise equal distances. That means R is not a Golomb ruler. 

Now assume that R is not a Golomb ruler. Thus, there are two equal 
differences |a — fe| = |c — d\ for some a, b,c,d£ R, where at least a, b and d are 
pairwise not equal. This means that {a, b, c,d} £ E and E ^ 0. D 

We can improve the running time of the construction to 0(|i?p) using a 
different approach. Instead of simply verifying every possible tuple, one can look 
at the distances between marks present in the ruler and examine which of them 
lead to edges in the graph. |HypergraphConstructionImproved| is a description of 
such an algorithm. 

In this algorithm we use an auxiliary map M to keep track of pairs of marks 
that measure specific distances. At first, we fill up this map: The first two loops 
iterate over distances present in R and add every pair of vertices to the entry in 
M corresponding to their distance. The map M then contains for every necessary 
distance in i? a list with all pairs of marks that measure this distance. In the 
second step, we add the edges to the designated characteristic hypergraph H: 
The last three nested loops again iterate over distances present in the ruler and 
simply add an edge to H for every pair of marks that measure this distance. 

For the correctness of this algorithm we first make some auxiliary observations: 
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Algorithm HypergraphConstructionlmproved: Constructing a clmr- 
acteristic liypergrapli for a given ruler 

Input: A ruler R cN. 

Output: A hypergraph H^ — {R,E). 

Start with an empty hypergraph H; 

Let R be the set of vertices in H; 

Create an empty map M that maps integers to lists; 

Smax <— max{a; : x € R} — min{a; : x e i?}; 

for i e i? do 

for jGR,i<j<i + *max/2 do 
I Add («, j) to the list mapped to j — i in M; 



8 for i G R do 



9 
10 

11 



for j ^ R,i < j < i + *moa;/2 do 

for {k, I) in the list mapped to j ~ i in M,j < k do 
^ Add the edge {i,j, k, 1} to H; 



12 return H: 



1 ^ 2^^^ 




'* 


*i 


1< i^ l" 

1 _ 2 ^rnax i 


1 



Figure 2.2: Shown is a ruler with four marks a, b, c, d. The pair a, c measures a 
distance that is above half the maximum distance S^^ax measurable by marks on 
the ruler. The pair b, d measures the same distance as a, c. However, because they 
measure such long distances, the measurements must overlap. By [Lemma 2.2.2"} 
there also is another distance (measurable without overlap) that also leads to 
this conflict. This distance must of course be below *max/2. 
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Lemma 2.2.2. Every edge in a characteristic hypergraph is due to two pairs of 
marks that measure the same distance and the measurements do not overlap. 

Proof. Assume that there are four marks a < b < c < d such that the following 
equation holds 

c— a = d — b=:S 

That is, the measurements of the pairs a, c and b, d overlap. Then subtracting the 
overlap c — 6 from the distance S gives b^ a = d — c, which implies measurements 
that do not overlap. D 

Corollary 2.2.1. Every edge in the characteristic hypergraph of the ruler R 
is due to a distance that is at most half the maximum distance measurable by 



marks on R (see Figure 2.2). 



Lemma 2.2.3. \Algorithm HypergraphConstructionImproved\ constructs a char- 
acteristic hypergraph for its input ruler. 

Proof. Assume that there are marks a<b<c<d in R such that {a, b, c, d} has 
not been added to H by the algorithm although the relation b—a = d—c< ^ma^/2 



holds. We can assume the properties of the marks because of Lemma 2.2.2 and 
[Corollary 2U\ 

It is clear that (a, b) and (c, d) have been added to M in lines 5 through 7. 
Then at some point in the execution of the algorithm, a must occur in line 8 
and b in line 9. Because (d, c) has been added to the list mapped to & — a in Af , 
the edge {a, 6, c, d} then is added to H. ^ D 



Lemma 2.2.4. 



Algorithm HypergraphConstructionlmproved runs in 0(|-R| ) time 



Proof. Obviously the running time is mainly dependent on the last three nested 
loops in lines 8-10. The two outer loops each iterate at most \R\ times. Retrieving 
the list from AI in line 10 can be done in 0(log(|i?|)) time, using red-black trees 
as implementation for M [3]. The iteration of the innermost loop in line 10 is 
bounded by a term in 0(|i?|), because any fixed distance S between two marks 
on the ruler R can occur at most 2|i?| times: S can be measured at most two 
times by one mark with any other mark. Adding the edge to the hypergraph 
is possible in 0(1) using an incidence graph representation of H and adjacency 
lists for example and thus the running time is in 0(|i?|^). D 

Note that we omitted long distances in the loop-headers in lines 6 and 9. 
However, the omission does not influence the asymptotic upper bound on the 
running time. This is a heuristic trick that could prove useful in practice. 

Lemma 2.2.5. The upper bound on the running time of \Algorithm Hypergraph] 
\ConstructionImproved\ is tight. 

Proof. There are characteristic hypergraphs that contain il(|i?p) edges: This 
holds for graphs constructed from rulers whose marks are taken consecutively 
from N. Consider the ruler R = {0, 1, ...,n} C N for even n; obviously every 
distance up to n is measured by marks in R. Let S < "■/2 be a fixed distance 
measured. How many possibilities are there to choose two pairs of marks such 
that they both measure S and the measurements do not overlap? One can place 
one pair leftmost on the ruler, count every possible placement of the other pair 
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to the right, then move the first one to the right by one and iterate. Summing 
over every distance up to "/2 this gives a lower bound on the number of edges 
in Hr: 

5=1 j=0 k=j+S 

No edge is counted twice here, because if there is an edge due to two different 
distances 1 < ^i < 62, the measurements of 62 overlap: let a < 6 and c < d both 
measure 61 with b < c. Where can S2 be placed? Certainly not between a, b and 
the other pair and also not between the pairs a, d and b, c. It can only be measured 
by both of the pairs a, c and 5, d and thus the measurements overlap. However, in 
the above construction, we are only counting non-overlapping measurements. D 

In conclusion, the following theorem now readily follows: 

Theorem 2.2.1. There is a hypergraph characterization for rulers such that 
Golomb rulers and only Golomb rulers correspond to hypergraphs without edges. 
The characteristic hypergraph for a ruler R can be computed in time 0(|i?p) and 
this bound is also tight. 

2.2.2 Notes on the Structure of the Characteristic Hyper- 
graphs 

At first, notice that the set of hypergraphs that can be constructed from rulers 



as in Subsection 2.2.1 ('Hypergraph Construction') is a strict subset of all 
hypergraphs with edges of size three and four. This is because the construction 
algorithm can be carried out using 0[n^) edge additions, n being the number of 
marks and thus vertices. However, generic 3,4-hypergraphs can have (?) G Q,{n'^) 
edges. 

2.2.2.1 Forbidden Substructures 

It would be very interesting to determine which sort of hypergraphs can and can- 
not be constructed; for example through a forbidden subgraph characterization. 
That is, a set F of hypergraphs, such that a 3,4-hypergraph if is a characteristic 
hypergraph for a ruler if and only if H does not contain a graph G G F in the 
sense of graph minors, subgraphs or induced subgraphs. Such a characterization 
would most likely be very useful in creating efficient algorithms for the problems 



we discuss in Section 2.4 (' Fixed- Parameter Tractability of Constructing Golomb 



Rulers ' ) . Unfortunately, we were not able to find a forbidden subgraph charac- 
terization. However, examples of forbidden and forbidden induced subgraphs 
have been found, which we discuss now. 



Lemma 2.2.6 (Small hand forbidden subgraph). The subgraph shown in Fig- 
\ure 2.3a\ cannot occur in a characteristic hypergraph. 

Proof. In an edge with three marks there is one mark exactly between the other 
two. Let a, b be two marks on a ruler. Where can a third mark in an edge already 
comprising b and a be? Either a, 6, or the new mark can be the mark in the 
middle, and thus there are at most three edges with three vertices intersecting 
in a and b. D 
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(a) Small hand (b) Large hand (c) Rotor 

Figure 2.3: Forbidden subgraphs 



Lemma 2.2.7 (Large hand forbidden subgraph). The subgraph shown in Fig 
\ure 2.3^ cannot occur in a characteristic hypergraph. 

Proof. In an edge with four marks there are two pairs of vertices that measure 
the same distance. We choose one unordered pair from {a, b, a} and thus define 
the distance that caused the edge. Then, for the fourth mark, there are only two 
possible positions left. Multiplying this with the number of possible unordered 
pairs, one gets six as an upper bound for such edges intersecting in three marks. 
However, every edge is counted twice here: Assume a < b has been chosen as 
pair. Then a fourth mark d is defined as follows: 

d = c-{b-a) V d = c + {b - a) 

These are the two possibilities for the fourth mark that are counted above. But 
we can reformulate these equations as follows: 

d = a+{c~b) V d = b+{c-a) 

These equations imply that the possibilities that are counted for the pair a, b are 
identified with one that is counted for the pair 6, c and one that is counted for 
a, c. This observation holds for all possible pairs from {a, b, c} and thus every 
possible location for d is counted twice. This means that there are at most three 
edges comprising four vertices that intersect in three vertices. D 



Lemma 2.2.8 (Rotor forbidden subgraph). The subgraph shown in Figure 2.3c 
cannot occur in a characteristic hypergraph. 

Proof. Fix a total ordering of the three marks in {a, c, d\. We then try to position 
b in that ordering and find that all possibilities lead to a contradiction. Because 
of the symmetry of the graph we can look at one specific ordering without loss 
of generality. 

Assume that a < c < d. Where can b be put? Assume that b < c. Because 
of the edge {6, c, d}, c is half-way between b and d. The edge {a, b, d} implies 
that either a = c (d) oy a < b. But then, because a, b are in one edge with c and 
in one with d, c — d. i The case 5 > c is symmetric. D 

Beside forbidden subgraphs, there are also some forbidden induced subgraphs: 
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(a) Scissors (b) Bird of prey 

Figure 2.4: Forbidden induced subgraphs 



Lemma 2.2.9 (Scissors forbidden induced subgraph). The graph shown in 
\Figure 2.4a\ is a forbidden induced subgraph in a characteristic hypergraph. 

Proof. We show that in such a configuration, an edge coinprising d, e and exactly 
one inark from {a, b, c} must also be present. 

We again use the fact that edges comprising four marks are due to two pairs 
of them having the same distances. Choose two pairs from {a, b, c} corresponding 
to the two edges. If those pairs comprise the same marks, the statement is trivial. 
If not, then the pairs must share one mark. Without loss of generality, let the 
pairs be a < 6 and b < c. The following equations hold; 



d = c±(6- 



e = a ± (c — &) 



Note that the sign before {b — a) cannot be negative at the same time with the 
sign of (c — b) being positive. Otherwise this would imply that d = e. In any 
other case, the two terms on the right hand side of the equations differ only in 
the sign of exactly two variables. This means that there exists an ?n e {a, 6, c} 
such that the following equation holds. 



m\ 



d\ 



Thus there is one edge {d, e, m}. 



n 



Lemma 2.2.10 (Bird of prey forbidden induced subgraph). The graph shown 
in \Figure 2.4b\ is a forbidden induced subgraph in a characteristic hypergraph. 

Proof. We show that there is a case distinction with two cases for 4-edges that 
contain two fix vertices a, b and if there are two 4-edges e^, Cj that contain a, b 
and correspond to the same case, then there is an edge CiU ej\{a,b}. 

This implies the statement of the lemma, because if there are three edges that 
intersect in two vertices then at least one additional edge must exist, making 
the graph in [Figure 2.4b| a forbidden induced subgraph. 

The case distinction is as follows: Fix a < & on a ruler. For the marks Ci and 
di of an edge e^ := {a, b, Ci, di} there are two cases: In the first case, both are 
between a and b or one to the left of a and one to the right of b. In the second 
case either both marks are positioned left of a or both marks are positioned 
right of b or one mark is positioned left of a or right of b and the other mark 
between the two. This distinction is exhaustive. 
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The existence of the additional edges can be seen as fohows: 
For edges e^ = {a,b,Ci,di} that correspond to the second case, it is clear 
that the marks in e^ \ {a.b} must have the same distance as a and b. Thus the 
existence of extra edges trivially follows in this case. 

Now assume the edges {a, 6, q, di} and {a, b, Cj, dj} are present in a charac- 
teristic hypergraph and correspond to the first case. Without loss of generality 
we assume Ci < di and Cj < dj. Then the following equations hold: 

a — Cj = dj — b a — Ci ~ di — b 

Subtracting the right from the left one, we get Ci — Cj — dj — di. This implies 
the edge {q, di, Cj,dj}. D 

These are some forbidden substructures we have found. There are other more 
specific sorts of forbidden structures, however, we could not observe a thorough 
explanation of the structure of the characteristic hypergraphs. We consider some 
general approaches to this question in the succeeding subsection. 

2.2.2.2 The Recognition of Characteristic Hypergraphs 

Abstracting from forbidden substructures one can pose the following question: 

Characteristic Hypergraph Recognition 

Input: A 3,4-hypergraph G. 

Question: Is there a ruler R such that Hr = G? 

We found two approaches to this question, which could lead to structural insights 
with proper analysis. We briefly describe them here. 

On the one hand, observe that any total order of the vertices in G defines 
for every edge or rather conflict a linear equation in the vertices or marks that 
are present in it. That is, any total order of vertices of G defines a homogeneous 
linear equation system with the vertices as variables, a solution to which can 
provide the sought ruler. However, certain solutions are infeasible: the trivial 
solution identifying every mark with zero, any other solution that requires two 
marks to be equal, and any solution that introduces edges that are not present in 
the graph. One challenge is to classify equation systems that exclude the second 
and third kind of solutions, another challenge is to either eliminate the need for 
a total ordering of the vertices or to find a way of deriving a subclass of feasible 
orderings in order for this classification to be efficient. A third challenge arises 
from the possibility of "derived edges" , that is, edges whose linear equations are 



linear combinations of other edges' equations (as we observed in Lemma 2.2.101. 
Another approach to Characteristic Hypergraph Recognition would 
be the following: Because every ruler i? is a subset of a ruler containing every in- 
teger between the maximum and minimum of R, one can investigate the question, 
whether there is an integer n such that G is isomorphic to an induced subgraph 



of i?{i, 2,. ..,«}• (See examples of such pseudo-complete graphs in Figure 2.5 ) 
This also would give a hint on whether G has a constructing ruler. 

An NP-hardness result for Characteristic Hypergraph Recognition 
of course would be a pessimal alternative. This would also mean that there 
is no classification through finitely many, constant-size, forbidden (induced) 
subgraphs. 
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(a) four vertices 



(b) five vortices 



Figure 2.5: Here, characteristic hypergraphs of the rulers {1,2,3,4} and 
{1,2,3,4,5} are shown. Since such graphs have the asymptotically maximum 
number of edges, we call them pseudo-complete. Any ruler of length < 5 has a 
characteristic hypergraph that is a subgraph of one of the graphs shown. Perhaps 
this property can be used to efficiently recognize characteristic hypergraphs. 




Figure 2.6: A subgraph of -ff{i.2,...,5} with three edges is shown. The incidence 
graph of this graph contains a K^^^ and thus is nonplanar. 



2.2.2.3 Other Properties 



Planarity and bounded degree are other properties that could possibly be ex- 
ploited for efficient algorithms. However, in general the characteristic hyper- 
graphs are neither: 

As we observed in [Subsection 2.2.1] there are characteristic hypergraphs 
with n vertices and J7(n^) edges. However, if a graph has bounded degree, i.e., 
there are at most c edges adjacent to any vertex, it has at most en edges for a 
constant c. 

A hypergraph is called planar if and only if its incidence graph representation 
is planar. We show an example of a hypergraph whose incidence graph comprises 
a complete bipartite graph with two vertex partitions of size three — a K^^^ — as 
a minor. Look at the three-edged subgraph of i?{i,2,...,5} shown in Figure 2.6 
All three edges intersect in 2, 3, and 4, and thus this graph's incidence graph has 
a -ftTa 3 as a subgraph. Obviously, this graph is then also part of the incidence 
graph of iJ{l,2,...,5}■ 
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2.3 A Simplified NP-Hardness Proof for Golomb 
Subruler 

Meyer and Papakonstantinou |i22l showed that the following problem is NP-hard 
via a reduction from a hard SAT- variant. 

Golomb Subruler 

Input: A finite set i? C N and I G N. 

Question: Is there a Golomb ruler R' C R with at least I marks? 

However, the construction of the ruler R corresponding to the SAT-formula 
is involved and hard to comprehend. We provide a reduction from a different 
SAT- variant, yielding a much simpler construction. In a simple corollary, we 
also gather a W[l] -hardness result for a similar problem, that did not directly 
follow from the original proof. To define the problem, which we reduce from, we 
first need the following definition: 

Definition 2.3.1. A boolean formula in 2-CNF is said to be antimonotone if 
and only if every literal is negative. 

Weighted Antimonotone 2-GNF SAT 

Input: An antimonotone boolean formula (p in 2-CNF and an integer k. 
Question: Is there a truth assignment for the variables in such that 
(f) is satisfied and at least k variables are assigned true? 

NP-Completeness. This problem is NP-complete via a straightforward re- 
duction from Independent Set. In Independent Set, a graph G — {V, E) 
and an integer k is given and it is asked whether there is a vertex set S* C 1^ of 
size at least k such that there is no e € i? with e C S. The set S is called an 
independent set of G. 

Let G = {V, E) be a graph. It is clear, that in any edge {u, v} £ E either u 
or V is not in an independent set of G. The following formula is true for any 
independent set S of G, where xs is the characteristic function of the set S. 

0(G) := f\ (-xs(u) V -xsM) 

{u,v}£E 

If we interpret the instances of x ^-s variables and 0(G) has a satisfying truth 
assignment with at least k positive variables, there is an independent set of size k 
for G and vice versa. 

Problem Restriction. We restrict the problem Weighted Antimonotone 
2-CNF SAT to instances such that (f) contains no two clauses consisting of 
the same variables. It is a simple observation that this has no impact on the 
complexity. (Note that such clauses cannot occur in the formula 0(G) in the 
construction above.) 

Reduction Outline. The basic idea of our reduction is to construct a ruler 
such that there are marks for every variable and additional marks that impose 
conflicts between those and only those marks that correspond to variables that 
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share a clause. In terms of characteristic hypergraphs, we construct a ruler R, 
such that Hfj contains only edges that consist of two marks that correspond 
to two variables and two auxiliary marks. Every edge shall correspond to one 
clause. 

The construction is given below. It is followed by some auxiliary statements 
that we use to show that the conflicts corresponding to the clauses are the only 
conflicts present in the constructed ruler. Once we have done that, we are able 
to prove that Weighted Antimonotone 2-CNF SAT can be solved with 

GOLOMB SUBRULER. 

Construction 2.3.1. Let <f) he an antimonotone boolean formula in 2-CNF with 
m clauses (numbered from 1 to m), let xi, ...,x„ he the variables in (p and k an 
integer. Construct a ruler R and an integer I as follows: 

Let V = {2('"+^)* : Xi in cj)} be the set of marks corresponding to the variables 
and let Ci = {2(™+2)* + 2^ — 1 : Xi occurs in the j 'th clause of (j)} and C = IJ^ Ci 
be the sets of marks corresponding to the clauses. 

The ruler R — VUC and I — k + 2m then constitute an instance of Golomb 

SUBRULER. 

Before we prove the polynomial running time and correctness, i.e., the 
equivalence of instances of this construction, we need some auxiliary statements: 

Lemma 2.3.1. The ruler /3^ := {2* : i e N} is a Golomb ruler. 

Proof. If /3^ is not a Golomb ruler, then there are integers i < j < k < I such 
that the following holds: 

2^' - 2^ == 2' - 2*= 

2J-i ^ 2*^-* - 1 = 2'"' 

However, because j, k < I, the right-hand side of this equation is strictly greater 
than the left-hand side, i D 

Corollary 2.3.1. The ruler V and each Ci is a Golomb ruler. 

Proof. The ruler F is a subset of /?^, and each Ci is a subruler of the following 
ruler: 

|2(m+2)z - 1 + j : j e ^N| p 

Now we know that the basic building blocks of the constructed ruler are 
conflict-free. We proceed to show that C also is conflict-free. Because conflicts 
do contain at most four marks, it suffices to prove that Ci U Cj U Cfe U C; contains 
no conflicts for any choice of i,j, k, I. We do that by successively adding more 
clause sets Co- 

While doing this we maintain that any one mark in V can be added to the 
unions of clause sets without imposing conflicts. Together with the observation 
that any one mark in a clause set can be added to V without creating a conflict, 
we then see that there are no conflicts in R that contain one mark from V and 
three of C or vice versa. 

To continue, we need the following auxiliary results that are a modification 
of lLemma 2.3.11 
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Lemma 2.3.2. For any integer d> 2 and any function / : N — > {0, ..., 2'^}, the 
following ruler is a Golomb ruler: 

/32(d,/):={2* + /(z):l<zGN} 

Proof. If there is such a ruler that is not a Golomb ruler, then there are four 
marks in it — i.e. i < j < k < I and constants < mi,mj,mk,mi < 2'^ — such 
that the following holds: 



(2* + m,) - (2* + m,) = {2'" + mi) ~ (2'"= + m^) 

If we divide this equation by 2^** , we get the following: 

(2''(J'-^) + TO,) - (1 + m,) = (2''('-*) + m) - (2'^('=-') + rrik) 
(2''(^'-^) + TOj) + (2'^('=-') + TOfe) - (1 + mi) = (2^('-*' + mi) 

Since < rhi, rfij, rfik, mi < 1 and j,k < I — 1, the left-hand side of the above 
equation is upper bounded by 2''('~*)~''+^ + 1. However, the right-hand side is 
lower bounded by 2'^^^~'^\ This is a contradiction, because the following relations 
hold: 1 < ; - i and 2 < d. i D 

Observation 2.3.1. The distance between two marks in (3^{d, f) is at least 2'^'^^^ . 

Proof. Let 2'^^ + mj and 2^** -I- m^ with z < j be two marks in /3J!(d, /). Then 
the following relations hold: 

2* -I- mj - (2* + mi) > 2'^i ~ [2'^' + 2'^) > 2''(*+i) - 2* - 2'^ 
= 2*(2''- 1) -2'' > 2''(2'' - 1) - 2'' = 2'^'^ -2"^+^ 

= 2 X 22"^-! - 2'*+! > 22''-i 

The last estimate holds because d > 2. D 

Lemma 2.3.3. The ruler 0^{d, f) does not cease to he a Golomb ruler, if we 
add the mark 2'^'^ + c to it, where i,c E N and < c < 2"^. 

Proof. It is clear that any edge e in the characteristic hypergraph of the 
ruler ^^(d, /)U{2*-f c} must comprise the marks a := 2* + /(z) and /3 := 2'^*-hc. 
The distance measured by those two marks is at most 2'^ and the distance mea- 



sured by any two marks in (3^{d, f) is at least 2 (Observation 2.3.1). Since 



d>2, this means that the edge e is not due to the distance |a — /3|. Because 



of Lemma 2.2.2 this also means that we can assume the marks a and /3 to be 
positioned between the other two marks in e. The existence of such an edge, 
however, has been disproven in [Lemma 2.3.2[ Have a look at the proof and 
observe, that we did not postulate the inequality of j and k and also did not 
make an assumption about the constants mj and m^. D 

Observation 2.3.2. Let v = 2(™+2)' eV, c^ e Ci and Cj,Cj e Q with i ^ j. 
The following relations hold: 

\v-c^\<2"' 2^"'+^ < \v - c,\ 

|c,-c,|<2™ 22™+^<|c.-c,| 
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Proof. The relations on the left are trivial. Those on the right-hand side follow 
from Observation 2.3.1 by modeling {v,Cj} and {ci,Cj} as subrulers of the 



ruler /3" (to + 2, /) with an appropriate function /. D 

Corollary 2.3.2. The rulers V -.^ V U {c}, c £ C and CJ := C-j U {v}, v eV 
are Golomb rulers for any integer j , 1 < j < n. 

Proof. The ruler V^ can be modeled as subruler of /3^(to, + 2, /) U {c} for 
appropriate values of /. According to [Lemma 2.3.3[ this is a Golomb ruler. 

Let V = 2'™+^)* and assume that C" contains a conflict. It is not possible 
that j = i, because otherwise the ruler Cj U {v} is a shifted subset of the Golomb 
ruler /3^. The case j y^ i is also not possible: A distanc e measured by c, , c , e Cj 
is at most 2™, however, it is 2^™+^ < jw — Cj| because of Observation 2.3.2 ^ D 



Lemma 2.3.4. The ruler C't U Cj is a Golomb ruler for all integers i,j with 
1 < i < j < n. 

Proof. Assume that there is an edge e in Hduc ■ This edge is due to two equal 
distances. A distance here can be measured by two marks in Ci (a), by two 
marks in Cj (b), or by one mark in Cj and one mark in Ci (c). 

Because of |Corollary 2.3.1] the edge e cannot be due to two distances both 



corresponding to (a), or (b). Because of Lemma 2.2.2 the case that both distances 



correspond to (c) reduces to the case that one distance corresponds to (a) and 
one to (b). 

Assume that the edge e is due to two distances corresponding to (a) and 
(b), respectively. Then there are k < I and o < p such that 2' — 2*^ = 2^ — 2° 



holds. We observed in Lemma 2.3.1 that /3 is a Golomb ruler. That means 
that / = p and k = o. Then, the variables Xi and Xj occur together in clause I 
and k. However, we excluded such instances from our considerations. 

Assume that the edge e is due to one distance di corresponding to (a) and 



one distance ^2 corresponding to (c). Because of Observation 2.3.2 the following 
relation holds. 

rfi < 2™ < 22™+3 < d2 

Thus, di and d2 cannot be equal, i 

The case (b) and (c) is analog to (a) and (c). D 

Corollary 2.3.3. The ruler {d U CjY := Cj U d U {v} is a Golomb ruler for 
all marks v ^V and integers i,j with \ < i < j < n. 



This follows directly from Corollary 2.3.2 and the proof of Lemma 2.3.4 



because every statement there still holds if we redefine Ci to contain v. 

Lemma 2.3.5. The ruler Ci U Cj U C/j is a Golomb ruler for all integers i,j,k 
with l<i<j<k<n. 



Proof. Assume that there is an edge e in HduC-uCk- Because of Lemma 2.3.4 
e must comprise marks from every of the three clause-sets. 

Assume that e is a 3-edge. It consists of one mark in C^, one in Cj and one 
in Cfe. However, these three marks form a ruler that is a subruler of /3^{d, f) 
with appropriate values for d and /. i 
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Now assume that e is a 4-edge. Because there are three clause-sets and e 
must comprise marks from every of them, e contains two marks mi and m2 from 
exactly one of the sets. Then, however, e forms a subset of a ruler corresponding 
to [Lemma 2 .3.31 i D 

Corollary 2.3.4. The ruler Ci U Cj U Ck U {v} with v E V is a Golomb ruler 
for all integers i,j, fc, I with l<i<j<k<n. 



This again directly follows from the proof of Lemma 2.3.5 

Lemma 2.3.6. The ruler C is a Golomb ruler. 

Proof. Because of [Lemma 2.3.4| and [Lemma 2.3.5] it only remains to show that 
Ci U Cj U Ck U C; is a Golomb ruler. This however is not hard to see, as any 
edge in the corresponding characteristic hypergraph must consist of four marks 
from every one of the clause-sets. Such an edge again would be a subruler 
of /3^{m + 2, /) with an appropriate function f.^ D 

Corollary 2.3.5. The ruler C U {v}, v £ V is a Golomb ruler. 

Corollary 2.3.6. // there are edges in Hr, there only are ^-edges that contain 
two marks from V and two from C. 

Lemma 2.3.7. // there is an edge e in Hr, then it intersects with V in two 
marks corresponding to the variables Xi and Xj and it intersects in exactly one 
mark with Ci and in exactly one mark with Cj . The variables Xi and Xj are 
together in one clause in <j). 

Proof. Let e be an edge in Hj^. Because of [Corollary 2.3.6| it intersects with V 
in two marks. Let these two marks be Vi = 2(™+^)* and Vj — 2('"+^)-', i < j, i.e. 
the marks corresponding to the variables Xi and Xj in (p. 

The edge e cannot contain a mark c^ € Ck with k ^ {i,j}, because then e 
would be a subset of a ruler corresponding to [Lemma 2.3.3| However, e can also 
not contain two marks Cfe,Cfc € Ck with k e {«, j}: Any distance measured by 
marks in {vk,Ck,Ck} is at most 2™ and those three marks are at least 2^™+"^ 
units away from the fourth mark (Observation 2.3^. 



Let e contain the mark Ci — 2(™+^^'' + 2'' — 1 e Ci and the fourth mark Cj £ Cj. 
Then, the following equations hold. 

c, - w, = Cj - Vj <^ Cj = 2("+2)J + 2''^ - 1 

This means that Xi and Xj are together in clause k. D 

Lemma 2.3.8. Let e, f be two different edges in Hji, then eD f C V. 

Proof. Assume that e, / have a non-empty intersection and e n / n Cfc 7^ 0. As 



we observed in Lemma 2.3.7 e n / also contains Vk £ V . This means that either 



e = / or there are three variables in one clause. ^ D 

We are now ready to prove the following: 

Lemma 2.3.9. Weighted Antimonotone 2-CNF SAT is polynomial-time 
many-to-one reducible to Golomb Subruler. 
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Proof. Let (p be an antimonotone boolean formula in 2-CNF with m clauses, 
let Xi, ...,x„ be the variables in (jj and let k be an integer, that is, let (/) and k 
constitute an instance of Weighted Antimonotone 2-CNF SAT. Construct 
an instance for Golomb Subruler, i.e., a ruler R — VUC and an integer k + 2m 
according to [Construction 2.3. 1| 

The marks to be constructed are given explicitly by the variables and their 
position in the clauses and the number of marks is polynomial in the number of 
variables and clauses. The length of the binary encoding of the marks is also 
bounded by a polynomial in the number of variables and clauses and thus the 
construction is possible in polynomial time. 

Given a truth assignment that satisfies </> with k positive variables, we can 
find a Golomb ruler with k + 2m marks in R: We have seen in [Lemma 2.3.71 that 
every edge in Hji corresponds to a clause. In one clause at most one variable can 
be assigned true. Therefore the marks corresponding to the positive variables 
form a Golomb ruler with k marks. Additionally, for every edge e in TJ/j, e n C 
contains only degree-one vertices, because of |Lemma 2.3.8| Because every edge 
contains at least one mark corresponding to a negative variable, we can add these 
degree-one vertices to the Golomb ruler constructed from the positive variables, 
yielding a Golomb ruler with k + 2m marks. 

Given a Golomb ruler R' C R with k + 2m marks, we construct a truth 
assignment that satisfies (j) with k positive variables: We assign every variable 
that corresponds to a mark in R' the value true, every other variable is assigned 
false. Because there are exactly n + 2m, marks in R, any such truth assignment 
has at least k positive variables. However, this simple assignment does not 
necessarily satisfy (j): There may be two marks in V and R' that are in the 
same clause in (j) and thus are in the same edge in Hji. But in every such 
edge e, there must be one mark r^ that is not in R' and, thus, we can simply 
exchange an arbitrary mark from e nV in R' with r^ and get an assignment 
that satisfies (j). D 

Now the following theorem directly follows: 

Theorem 2.3.1. Golomb Subruler is NP-complete, even if there are only 
conflicts with four marks in the input instance. 

The reduction also yields a W[l]-hardness result for a modified problem: 

Golomb Subruler > Double Conflicts 

Input: A ruler i? C N and fc e N. 

Question: Is there a Golomb ruler R' C i?, such that \R'\ is at least k 

plus two times the number of edges in Hal 



Corollary 2.3.7. Golomb Subruler > Double Conflicts is W[l]-hard 
with respect to parameter k. 

Proof. Observe that there is a parameterized reduction from Independent 
Set parameterized with the sought independent set size to Weighted Anti- 
monotone 2-CNF SAT parameterized with the number of positive variables 



in a satisfying truth assignment (we have given the idea for this on page 21 ). 
The reduction from Weighted Antimonotone 2-CNF SAT to Golomb 
Subruler we have given above maps a formula (p and the parameter "number of 
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positive variables" fc to a ruler R and the sought ruler size k + 2m, where m is the 
number of edges in Hn. This means, this reduction identifies the parameters of 
Weighted Antimonotone 2-CNF SAT and Golomb Subruler > Double 
Conflicts, making it a parameterized reduction. D 



2.4 Fixed- Parameter Tractability of Construct- 
ing Golomb Rulers 

The number of implementations of the search for optimal or near-optimal Golomb 
rulers hints to the importance of this problem (see e.g. [TU l [TT | [T^ [T ^ [251 [2^ [55] ) . 
Unfortunately, several problems closely related to the construction of Golomb 



rulers have been proven to be NP-complete (cf. Subsection 1.3.3 ('Complexity 



Theory] ')). 

We now look at the fixed-parameter tractability of constructing Golomb 
rulers with the goal to get exact and relatively efficient algorithms. In particular, 
we focus on the following problem. 

Golomb Subruler 

Input: A finite set i? C N and n e N. 

Question: Does there exist a Golomb ruler R' C R oi at least n marks? 

For a given R in Golomb Subruler, let Hr = {R, E) be the characteristic 
hypergraph as defined in Subsection 2.2.1 ( '[Hypergraph Construction ' ) . Now 



the task is to find a subset R' C R such that the induced subgraph Hji> = Hji[R'] 
contains no edges. This can either be done by selecting marks to keep or by 
deleting marks from R. 



2.4.1 Mark Deletion Parameter 

If we decide to delete marks, we can parameterize Golomb Subruler with 
the maximum number of allowed mark deletions. Together with the notion of 
hypergraph characterization, we can reformulate it as follows. 

Golomb Subruler Mark Deletion 

Input: A finite set i? C N and fc e N. 

Question: Is there a set of marks D with size at most k, such that 

Hii\D contains no edges? 

The above described graph problem is strikingly similar to the Hitting Set 
problem. In HiTTiNG Set, a hypergraph is given and a (minimum size) subset S 
of vertices is sought, such that every edge has at least one vertex in S. In fact, 
our problem can canonically (by simply computing the characteristic hypergraph) 
be reduced to this problem parameterized with the size of S in polynomial time. 
However, as we observed in [Subsection 2.2.2] our graph instances are a 
strict subset of all possible hypergraph instances, which raises hope for better 
algorithms than those known for this generic problem. 
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Algorithm SearchTreeNode: Solving Golomb Subruler Mark Dele- 

TION. 

Input: A liypergrapli H = {R, E), an integer k and a vertex deletion 

set D. 
Output: A set of vertices D of size < k such that H[R\ D] contains no 

edges, if it exists. 
if H[R\ D] has no edges then Output D and halt; 
if \T\ = k then Abort this branch; 
else 

Choose an edge e £ E; 

for i G e do SearchTreeNode (H, k ~1,DU {i}) 

Algorithm SearchTree: Solving Golomb Subruler Mark Deletion. 
Input: A ruler R cN and an integer k. 
Output: A set of marks D of size < k such that R\ D is Golomb, if it 

exists. 
Hji ^— HypergraphConstructionCi?) ; 
SearchTreeNode (i/i{, fc, 0) 



2.4.1.1 Fixed-Parameter Algorithm 

From the hypergraph characterization and the notion of deleting vertices we 
immediately get a search tree algorithm: It is clear that in every edge of a 
characteristic hypergraph at least one mark has to be deleted. That means, one 
can choose one edge, branch into all possibilities of deleting one mark in that edge 
and do this recursively until either k marks have been deleted or the characteristic 



graph has no edges. (See also the pseudocode in Algorithm SearchTree ) 



Because every edge in a graph defined by a ruler R as in Subsection 2.2.1 



( '[Hypergraph Construction ' ) has edges with at most four vertices, and the 



recursion depth of SearchTreeNode| is at most fc, the running time of this 



algorithm is time-bounded by a term in 0*(4'^). That is, this algorithm suffices 
to classify this problem as fixed-parameter tractable with respect to at most k 
mark deletions. 

As we noted above, Golomb Subruler Mark Deletion can be solved with 
Hitting Set algorithms. The fastest known algorithm for Hitting Set with 
edges of four vertices and parameterized with the solution size k is 0(3.076'^ -|-m) 
by Dom et al. [15]. For this more restricted problem it is probably possible 
to beat this upper bound. Unfortunately we did not find such an algorithm. 
However, there are some applicable heuristic tricks that improve the above 
mentioned search strategy in practice: 

Edges with Three Vertices. In every edge with three vertices, at least one 



vertex has to be deleted, to make the graph edge-free. Thus, in SearchTreeNode 
one can search for such an edge, and branch into the deletion of every one vertex 
in it. This leads to a branching vector of (1, 1, 1) and the search for the edge 
takes 0(m) time, m being the number of edges in the graph. However, there are 
characteristic hypergraphs that do not comprise edges with three vertices and 
thus this rule does not suffice to improve the theoretical upper bound for the 
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search running time. 

Dominating Vertices. Another simple rule for an improved branching vector 
is the search for dominating vertices. 

Definition 2.4.1 (Dominating vertices). A vertex v dominates a vertex u, if 
and only if v is in every edge that u is incident to. 

It is clear, that if there is an optimal solution that contains a vertex u that 
is dominated by another vertex w, there is also an optimal solution that contains 
V instead of u (and obviously an optimal solution does never comprise both). 



Thus, in SearchTreeNode one can search for a vertex v that dominates another 
vertex, choose an edge that contains both and branch on either deleting v or any 
one of the vertices in this edge, that are not dominated by v. The worst case for 
this rule is, that there are only vertices that dominate exactly one other vertex, 
leading to a branching vector of (1, 1, 1). The search for dominating vertices can 
be conducted in time 0{nm?) by iterating over every vertex v and checking for 
every adjacent vertex, whether its incident edge set is a subset of the edge set of 

V. 

Also notice that this rule doubles as a reduction rule for edges that intersect 
in at most one fix vertex v with any other edge in the graph. 

Cementating Branched-on Vertices. A strategy that is applicable to any 
deletion search is the "cementating" of an entity w, when the recursive call for 
the deleted v returns and no solution has been found. It is then clear, that no 
solution with v deleted is possible in further branching and v can be prohibited 
from deletion or be "cementated" . 

This strategy however can be extraordinarily powerful in the search routine 
for our particular problem: Suppose that at a call of |SearchTreeNode| a set 
of vertices C has been cementated. This means, that if there is a solution to 
be found through further branching, this solution does not contain any vertex 
in C, and C is a Golomb ruler. This can already be exploited for an early 
termination of the branch, if C is not. Furthermore, it is clear that any distance 
between vertices in C must not appear in the rest of the vertices. This has 
two implications: First, if a distance appearing between vertices in C is also 
measured by one vertex w in C and one vertex u not in C, then u has to be 
added to the solution. Second, any such distance measured by two vertices not 
in C has to be destroyed, leading to a branching vector of (1, 1). 



This rule is very powerful in practice as we observe in chapter 3 ( ' Implemen- 



tation and Empirical Results ' ) and with further analysis might also lead to an 



improvement of the theoretical upper bound on the running time of the search. 

2.4.1.2 Cubic Problem Kernel 

Abu-Khzam [Ij observed, that Hitting Set, parameterized with the maximum 
number of vertices d in an edge and the maximum solution size fc, has a problem 
kernel of at most 0{¥^~^) edges. However, the reduction rules are not directly 
applicable to our problem. This is because instances produced by those reduction 
rules may not correspond to a ruler anymore: Some rules shrink edges, but 
since the edges correspond to conflicts in characteristic hypergraphs, there is no 
equivalent to a shrunk edge. 
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However, there are adequate substitute rules for our hypergraph instances that 
also retain the problem kernel size of 0{k^) edges (in characteristic hypergraphs 
d — 4). These substitute rules compared to the ones by Abu-Khzam [1], we 
salvage the basic idea of the high occurrence rules and use some structure we 
observed in characteristic hypergraphs to make them applicable. 

In the following we assume that the characteristic hypergraph of the input 
ruler has been computed, and is kept updated alongside the ruler. We first 
list some simple and obvious rules. Then we give a rule that suffices to bound 
the number of 3-edges in the characteristic hypergraph. With an additional 
observation we give another rule to bound the number of 4-edges in the reduced 
graph. With the help of these two bounds, we are then able to bound the number 
of marks in a reduced instance. 

Reduction Rule 2.4.1 (Lone edges). If there is an edge that does not intersect 
with any other edge, then remove it and all vertices it comprises from the graph 
and reduce k by one. 

Reduction Rule 2.4.2 (Lone vertices). // there is a vertex with degree zero, 
then remove it from the graph. 

Reduction Rule 2.4.3 (Leaf edges). // there is an edge that intersects any 
other edge in at most one of its vertices v, then remove v, remove any edges 
incident to v and decrement k. 

It is clear, that these three simple rules are correct and can be executed in 
time 0(n + m). 

Reduction Rule 2.4.4 (High Degree for 3-Edges). // there is a vertex that 
has more than 3k incident 3-edges, then remove it from the graph, remove any 
incident edges and reduce k by one. 



Lemma 2.4.1. Reduction Rule 2.4.4 ** correct and can be carried out in running 
time 0{n + m). A graph has at most 3fc^ 3-edges, if this rule cannot be applied 
to it and it is solvable with k vertex deletions. 

Proof. Assume there are more than 3fc 3-edges incident to one vertex v. We 
have seen in [Subsection 2.2.2[ [Lemma 2.2.6| that there are at most three edges 
with three vertices that intersect in two vertices. That means the deletion of any 
other vertex in the graph can destroy at most three edges incident to v. Thus, if 
V is not deleted, at least fc + 1 vertices are necessary to destroy all edges incident 
to V. 



To apply Reduction Rule 2.4.4 one can simply iterate over every 3-edge and 
count the occurrence of the vertices in an array indexed by the vertices. This is 
possible in time 0(n + m). 

Now assume pleduction Rule 2. 4. 4| cannot be applied to the yes-instance H. 
Every edge in H has to be destroyed, 3-edges in particular. One vertex can hit 
at most 3k 3-edges and thus, H has at most 3k^ 3-edges. D 

Lemma 2.4.2 (Induced Clique Structures). If there are more than 3k 4-edges 
that intersect in two vertices, this instance cannot be solved with at most k vertex 
deletions. 
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Proof. Recall that if there are more than two 4-eclges intersecting in two vertices, 
there are additional edges not containing those vertices, as we observed in 



Lemma 2.2.10 in Subsection 2.2.2 In fact, in this lemma we proved a stronger 
result: There is a case distinction with two cases for 4-edges that intersect in 
two vertices a, b and if the edges ei and 62 correspond to the same case, there is 
an edge (ei U 62) \ {a, b}. 

In case one, the marks other than a, 6 are in between a and b or one to the 
left of a and one to the right of b. That means, edges that correspond to case 
one cannot intersect in vertices other than a, b and if there are more than fc — 1 
such edges, there is a structure with more than k pairs of vertices that are all 
pairwise in an edge. It is clear that in such a structure, more than k vertices 
have to be deleted, to make it edge-free. 

In case two, the two marks other than a, b are positioned left of a or both 
marks are positioned right of b or one mark is positioned left of a or right of b 
and the other mark between the two. In this case a, b measure the same distance 
as the other two marks in the edge. That means, at most two edges in this case 
can intersect in a vertex other than a and b. Thus, if there are more than 2fc 
edges that intersect in two vertices and correspond to case two, then there is a 
structure with more than 2fc + 1 pairs of vertices that are all pairwise in an edge. 
At most two of those pairs overlap in one vertex, and thus more than k vertices 
have to be deleted to destroy every edge in this structure. 

If there are more than 3/c edges that intersect in two vertices, either more 
than 2k edges correspond case one or more than k to case two and thus the 
graph cannot be solved with k vertex deletions. D 

Reduction Rule 2.4.5 (High Degree for 4-Edges). // there is a vertex that has 
more than 3k^ incident 4-edges, remove it from the graph, remove any incident 
edges and reduce k by one. 

Lemma 2.4.3. \Reduction Rule 2.4. 5\ is correct and can be carried out in time 
0(n + m). A graph has at most 3fc'^ 4-edges, if this rule cannot be applied to it 
and it is solvable with k vertex deletions. 

The proof is analogous to the proof for |Lcmma 2.4.1| by substituting |Lcmma 2.2.6| 
with lLemma 2.4.21 

Theorem 2.4.1 (Problem Kernel for GOLOMB Subruler Mark Deletion). 
GOLOMB Subruler Mark Deletion has a problem kernel with at most 
9k^ + 2fc^ + k marks. The characteristic hypergraph of the ruler of a kernelized 
instance has at most 3fc^ 4-edges and 3fc^ 3- edges. The kernel can be computed 
in 0{kn + km) time if the characteristic hypergraph is known. 



Proof. To compute the kernel, proceed as follows: Apply [Reduction Rule 2.44 



apply [Reduction Rule 2.4.5| and recurse until neither applies anymore. Then 



apply [Reduction Rule 2.4T2 until it does not apply anymore. 



Since both high-degree reduction rules can be applied at most k times, the 
procedure recurs at most k times and the overall running time adds up to 
0{kn + km). [Reduction Rule 2.4.2 can of course be applied exhaustively in 



0(n) time. 

The upper bound on the 3- and 4-edges follows from [Lemma 2.4. 1[ and 
[Lemma 2.4.3[ In a yes-instance there is a set S of at most k vertices, such that 
every edge in the graph has a non-empty intersection with S. That means in 
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each of the 4-edges, there are at most three vertices not in S and thus there are 
at most 9k^ + k vertices in 4-edges. This argument holds analogously for 3-edges 
and thus there are at most 9k^ + 2k^ + k vertices in a yes-instance. D 

2.4.2 Mark Preservation Parameter 

We also tried to analyze Golomb Subruler with respect to the size of the 
sought Golomb ruler. However, the problem mostly escaped our attempts. We 
conjecture it to be W[l]-hard, and gather some hints towards this in the following. 
With the notion of hypergraph characterization, Golomb Subruler refor- 
mulates as follows. 

Golomb Subruler 

Input: A ruler i? C N and fc e N. 

Question: Is there a ruler R' C R with size at least k, such that Hn/ 

contains no edges? 

At first, observe that GOLOMB Subruler can be solved with Independent Set 
on hypergraphs. In INDEPENDENT Set on hypergraphs, a hypergraph H is given 
and a (maximum size) vertex set S is sought, such that H[S] contains no edges. 
However, Independent Set on r-uniform hypergraphs, parameterized with 
the size of the sought vertex set, has been proven to be W[l]-hard by Nicolas 
and Rivals [25j . Their proof relies on heavily overlapping edges, and thus this 
approach is not directly applicable to our problem. Nevertheless, our problem 
retains some of the features of general Independent Set. For example, a 
simple branching strategy like the one in Hitting Set seems not to be feasible, 
because in our instances too, there can be edges that do not contain solution 
vertices at all. 

If we modify the problem slightly, it indeed becomes W[l]-hard: 

Golomb Subruler > Double Conflicts 

Input: A ruler RCN and fc e N. 

Question: Is there a ruler R' C i?, such that Hji' contains no edges 

and i?' has a number of marks that is least k plus two times the number 

of edges in Hj^l 

The proof for this is given in [Section 2.3| on page [26] 
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Chapter 3 

Implementation and 
Empirical Results 

In this chapter, we investigate the practical imphcations of our considerations in 



Section 2.4 ('Fixed-Parameter Tractabihty of Constructing Golomb Rulers'). 



We implemented the search tree algorithm that is discussed there. It uses some 
heuristic improvements and the problem kernel we observed. We first describe 
the implementation in detail and then report on our results. Some of the heuristic 
improvements prove very effective in practice. 

3.1 Description of the Implementation 

Problem Definition. Our implementation solves the following problem with 
help of the fixed-parameter algorithm we developed in [Subsection 2.4. l[ 

Maximum Mark Golomb Subruler 

Instance: A ruler R cN. 

Task: Find a Golomb ruler _R' C i? of maximum cardinality. 

Observe that, theoretically, this problem can be used to answer a variety of 
questions: By using appropriate rulers — with marks taken consecutively from N — 
the following problem reduces to the one above. 

Instance: An integer D E N. 

Task: Find a Golomb ruler of length at most D and maximum number 

of marks. 

Golomb Ruler Decision can be solved with the second problem, which 
means that our algorithm can be used to directly search for optimal Golomb 
rulers as well as a subroutine in other programs that produce partial instances 
corresponding to Maximum Mark Golomb Subruler. 

Algorithm. The algorithm first generates the characteristic hypergraph of the 
input ruler R. It then proceeds to heuristically compute a minimal subset S of 
R such that Hj^^g contains no edges. This solution is used as an upper bound 
for the parameter in the corresponding instance of Golomb Subruler Mark 
Deletion. We interpret this upper bound and the input ruler as an instance 
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of GOLOMB SuBRULER MARK DELETION, solve it, and Successively decrease 
the parameter, thus obtaining successive smaller mark deletion sets, until no 
solution can be found anymore. The last solution S that could be found is a 
minimum-size set of marks such that Hji contains no edges. It is clear that 
R\ S must be optimal to Maximum Mark Golomb Subruler. See also 
|Algorithm FindGolombSubruler] 

Algorithm FindGolombSubruler: Computing Golomb subruler with 
maximum number of marks for a given ruler 
Input: A ruler i? C N. 

Output: A Golomb ruler R' Q R with maximum number of marks. 
Hn <— HypergraphConstructionlmproved(iZ); 
MinimalMarkDeletionSet -(— SolveHeuristic(iJij); 
continue -(— true; 
A; ^ |MinimalMarkDeletionSet|; 

AtMostfcMarkDeletionSet ^ MinimalMarkDeletionSet; 
while continue do 

ParameterizedSolution ^ SolveParcuneterizedCff/j^, k. 0, 0); 
if ParameterizedSolution = false then continue <— false; 
else 

k -i— |ParameterizedSolution| — 1; 
AtMostfcMarkDeletionSet -(— ParameterizedSolution; 

R' ^ R\ AtMostfcMarkDeletionSet; 
return R'\ 



The heuristic solution is computed via two very simple strategies and the 
best solution is kept. The first strategy greedily deletes vertices from the input 
graph Hr, until it is edge- free. The second strategy greedily selects an edge, 
deletes all vertices of this edge, and iterates until the graph contains no edges 
anymore. Observe that the second strategy yields a solution that contains at 
most four times the number of vertices in an optimal solution, because in any 
edge at least one vertex has to be deleted. 

For the subroutine ISolveParameterizedl that solves the instances of Golomb 
Subruler Mark Deletion, we make use of the (heuristic) improvements 



described in Subsection 2.4.1.1 That is, we cementate vertices (see page 291, if 
the recursive calls of |SolveParameterized| return without finding a solution, we 
use the implications of the cementated vertices and we try to find edges that 
imply favorable branching vectors. Additionally, in search tree nodes, we apply 
the reduction rules we described in|Subsection 2.4.1.2 A pseudocode description 



of the search tree subroutine can be seen in Algorithm SolveParameterized (The 



cementating of marks in lines 10 through 14 is simplified for readability.) 

Implementation. We implemented the above described algorithm in the 

Objective Caml language. Objective Caml is a multi-paradigm language, allowing 
for object oriented, functional and procedural programming styles. We chose it 
because first, it allows for easy transition from theoretical algorithms to practical 
programs via the functional programming possibilities. Second, because it is 
very robust to programming errors through strict variable typing and automatic 
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Algorithm SolveParameterized: Solving Golomb Subruler Mark 

Deletion 

Input: A characteristic Hypergraph Hji, an integer fc, a set of cementated 

marks C C R and a solution set S. 
Output: A set of marks S such that \S\ < k and Hp\g contains no edges 
or false if such a set does not exist. 

1 if fc < then return false; 

2 if Hji contains no edges then return S; 

3 Apply the problem kernel to Hji and fc, adding deleted vertices to S; 

4 if fc < or Hii exceeds the problem kernel size then return false; 

5 if Hr contains no edges then return S; 

6 BranchingEdge -s— An edge e in Hj^ such that e\C has minimum size; 

7 for i S BranchingEdge \ C do 



8 
9 
10 
11 
12 
13 
14 



Solution -(— SolveParcmieterized(_ffjj.\ r-ji , fc — 1, C, SU{i}); 
if Solution = false then 

C^CU{i}; 

for a, b,c E C do 

X ^ {c+\a-b\,c-\a-b\}; 

if A n C 7^ then return false; 

Hr^ Hn\x\ S ^ S U (X n R); k ^ k ~ \X n R\; 

else return Solution; 



16 return false; 



type inference. And third, because the implementation is meant to be a proof of 
concept and not a highly optimized solver. 

The design is a simple one-process, one-threaded solution. 

Testing Environment. The experiments were conducted on an Intel Xeon 
E5410 machine with 2.33 GHz, 6 MB L2 cache, and 32 GB main memory. The 
operating system was GNU/Linux with a kernel of version 2.6.26. The Objective 
Caml compiler used was of version 3.10.2. We did not utilize the symmetric 
multiprocessing capabilities available on the system. 

3.2 Running Times and other Results 

We conducted our experiments on rulers with n marks of the form {i G N : i < n}. 

General Running Times. Varying the number of marks n, we observe the 
running time behavior shown in [Figure 3.1| The running times obviously 
increase exponentially with increasing n. However, if we fit the running time 
function f{n) = a + bri^c" to the empirical running times via the Levenberg- 
Marquardt algorithm [24], we get the following result: 

f{n) = 0.989502 + 9.89765 x IO^^ti^I. 18727" 

Even with iteratively using the 0*(4'^) subroutine, we still obtain an algorithm 
whose empirical running time has an exponential term below 1.2". Although 
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Figure 3.1: Observed running times of Algorithm FindGolombSubruler 
number of marks of the input ruler. 
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Figure 3.2: Greedy solution size minus optimal mark deletion set size, 
greedy algorithm performs very close to the optimal solution sizes. 



The 



this seems like a good result, it would be necessary to examine rulers of about 
550 marks, to proof or disproof the optimality of the currently shortest known 
Golomb ruler with 27 marks (as of now, Golomb rulers of minimum length with 
up to 26 marks are known). Assuming the running time behavior corresponds to 
the function / we determined above, this would result in running times in the 
order of lO'^" years on our architecture, i.e., this particular implementation is 
not feasible for this application. However, with the reduction rules and heuristic 
improvements, it still might be possible to use it as subroutine in such search 
algorithms. Experiments on random rulers as input could settle this question. 



Greedy Strategies. In Algorithm FindGolombSubruler two greedy strategies 
are used to obtain an upper bound for the parameter. We observe that in our 
experiments the strategy that greedily deletes vertices always yields a solution 
that is superior to greedily deleting all vertices of one edge. The size of the 
greedy solution is very close to the optimal solution in our experiments, as shown 
In this context, it would be interesting to test whether there are 



Figure 3.2 



lower bounding techniques that can be applied fast. 



Cementating Vertices. The cementating of vertices that have a recursive 



call returning without solution yields large speedups as is shown in Table 3.1 



36 



2,734 


3,278 


893,070 


1,351,977 


60,035,055 


103,541,222 



Number of Marks No Cementation Cementation 

20 35.22 0.18 

25 8313.78 1.27 

Table 3.1: Cementating vertices yields tremendous speedups. Running times in 
seconds. 



Number of Marks Search Tree Nodes Vertices Deleted 

20 
40 
60 

Table 3.2: Search tree nodes and vertices deleted due to cementation in compari- 
son. 



In [Table 3.2[ it is shown, how many vertices are deleted by it in the course of 
the search tree, when varying the size of the input ruler. This success can be 
explained with the following example: 

Suppose branching has been done on edges that contain vertices of a confined 
region on the input ruler. A set of d cementated vertices have been gathered and 
now the algorithm moves on to branch on an edge in a different region on the 
ruler. If the first recursive call for a vertex in the edge returns negative, d{d— 1)/2 
vertices can be deleted, because the vertex now is cementated and the set of 
cementated vertices must be a Golomb ruler. The number of deleted vertices 
most likely increases another time, when the recursive call for the second vertex 
in the edge returns negative and so on. This means that we get a branching 
vector of (1, d{d — l)/2, d{d — l)/2 + c, ...) for this particular case. The running 
times suggest that the branching vector of an edge is likely to be much more 
favorable in practise with cementated vertices. 

Reduction Rules in the Search Tree. When not using cementation, we 
obtain speedups of about two using the high-degree rules. However, when 
using cementation they do not yield running time benefits anymore; although 
the reduction rules delete two to three vertices per search tree node and are 



responsible for many direct terminations of branches (as we show in Table 3.3). 
The time needed to calculate the rules outweighs their benefits in our instances. 
This is shown for the high-degree rule in [Table 3.4[ This can not be remedied by 
only applying the reduction rules every d'th search tree node in our instances. 
However, considering the number of deleted vertices and aborted branches and 
the fact that they seem to improve in overall effectiveness with increasing number 



of vertices (see Table 3.3), it might prove worthwhile to thoroughly optimize the 



implementation to support these routines in minimal time. 

Conclusion On the plus side, we could show that our heuristic improvements 
are very effective and our reduction rules delete many vertices and are effective 
in trimming the search tree. Unfortunately, despite a relatively small exponential 
term in the observed running times, this did not lead to an implementation that 
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Number 

of Marks Search Tree Nodes Vertices Deleted Branches Aborted 

20 2,734 6,148 1,856 

40 893,070 2,472,270 570,670 

60 60,035,055 186,900,562 34,506,053 

Table 3.3: Search tree nodes compared to vertices deleted and branches aborted 
due to the high-degree, leaf-edge and lone-edge reduction rules. 



Number 


With 


Without 


of Marks 


High-Degree Rules 


High-Degree Rules 


20 


0.18 


0.17 


30 


6.07 


5.48 


40 


106.21 


97.79 


50 


1212.30 


1134.45 



Table 3.4: Running times using and ignoring the high-degree reduction rules 
whilst cementating vertices. All values in seconds. 



is feasible for discovering new Golomb rulers. 
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Chapter 4 

Conclusion and Outlook 



Our Work. In this work, we have given a short overview of some hurdles 
and algorithmic possibilities in the field of Golomb rulers. The basic problems 
Optimal Golomb Ruler and Golomb Ruler Decision seem to be elusive 
to classic complexity classification and may well lie between P and NP. To settle 
these questions would of course be very interesting; however, it seems likely 
that this would imply major new insights into classic complexity theory and/or 
number theory. 

The natural hypergraph characterization we have given for rulers makes it 
possible to get a better understanding of conflicts with respect to Golomb rulers. 
We observed some structure in characteristic hypergraphs that we later exploited 
for a problem kernel. This implies that a more sophisticated structurization 
of those graphs could lead to much better algorithms and thus would be very 
interesting. A more thorough understanding of the graphs could also be used 
to settle other questions related to Golomb ruler construction, for example the 
W[l]-hardness of Golomb Subruler with respect to the mark preservation 
parameter. 

We have given a simplified proof for the NP-hardness of Golomb Subruler 
that also lead to a W[l]-hardness result for a modified problem that did not 
directly follow from the original proof. We hope that this makes fixed-parameter 
research for related problems more accessible and attractive. 

Concerning fixed-parameter algorithmics, we have observed that Golomb 
Subruler Mark Deletion is tractable and we have given a cubic-size problem 
kernel. We strongly believe that this bound can be surpassed and we also have 
just scratched the surface regarding solution algorithms. 

The implementation of a corresponding algorithm showed that the reduction 
rules prove quantitatively efficient in a search tree, however, their effective 
running times have to be improved. The technique of cemcntating vertices 
applied to Golomb rulers proved very effective in practice. 

Other Aspects. There also are some interesting aspects which we did not 
cover in this work. We briefly list some of them here: 

• To our knowledge the relations between Difference Cover, Turnpike 
and Optimal Golomb Ruler have mostly been left unexplored in the past. 
There might be a trinity similar to Clique, Vertex Cover and Independent 
Set. 
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• Also to our knowledge the parameterized complexity of GOLOMB Ruler 
Decision has not been touched yet. Since easily checkable lower bounds and 
allegedly efBciently constructible upper bounds on the length of a Golomb ruler 
exist |12j . it may be possible to apply a below or above guarantee parameteriza- 
tion. 

• There is a problem called Minimum Redundancy Linear Array |23j , 
where rulers are sought which are perfect — they measure every distance up to 
their length — and have minimum redundant differences between marks. Since 
in a Golomb ruler there is no redundancy at all, maybe this would make for a 
practical parameterization. 

• Considering the research in local search and evolutionary strategies (TUl HI] 
[551 1311 133] , it might be interesting to contemplate fixed-parameter tractability 
for problems arising in this field. 

• The notion of treewidth is a popular topic in fixed-parameter algorithms. 
There is a theorem by Courcelle [18] that classifies every decision problem that can 
be formulated in monadic second-order logic as fixed-parameter tractable with 
respect to the parameter treewidth. One could investigate if such formulations 
exist for the hypergraph problems we considered in this work. 
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